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USING ROBOTIC FISH TO EXPLORE THE HYDRODYNAMIC MECHANISM 
OF ENERGY SAVING IN A FISH SCHOOL 

LIANG LI, LICHAO JIA, AND GUANGMING XIE 


Abstract. Fish often travel in highly organized schools. One of the most quoted func¬ 
tions of these configurations is energy savings. Here, we verified the hypothesis and ex¬ 
plored the mechanism through series of experiments on “schooling” robotic fish, which can 
undulate actively with flexible body, resembling real fish. We find that, when the school 
swims in the same spatial arrays as the real one, the energy consumption of the follower 
mainly depends on the phase difference, a phase angle by which the body wave of the 
follower leads or lags that of the leader, instead of spatial arrays. Further analysis through 
flow visualization indicates that the follower saves energy when the phase difference con'e- 
sponds to the situation that the follower flaps in the same direction of the flow field induced 
by the vortex dipole shedding by the leader. Using biomimetic robots to verify the biolog¬ 
ical hypothesis in this paper also sheds new light on the connections among the fields of 
engineering, physics and biology. 


1. Introduction 

In the animal kingdom, many species move in highly organized formations IHEIEI. 
Some well-known systems include the queuing of migrating lobsters , organized swarm¬ 
ing ducklings 0, flocking birds 0, and schooling fishes Ena. Apart from the social 
benefits Ellllll and physiological factors im, energy saving for individuals is another 
prevalent hypothesis of why animals prefer to move in groups. In the case of fish schools, 
since Weihs first predicted that fish swimming in school will save energy d, this hypoth¬ 
esis has been proved through series of methods, such as theoretical analysis ini, compu¬ 
tational methods m Ha mini, biological experiments experiments mmniioi® and 
physical models ETl l22l |2^ . For theoretical analysis and computational methods, turbu¬ 
lent flow is difficult to estimate. These methods can only draw an approximate conclusion. 
For experiments on real fish, it is difficult to control the fish school and decouple the social 
factors from physical factors. Recently, with the development of mechanics, electronics, 
materials and controls, it is possible to build a man-made ’’fish school” to unravel the hy¬ 
drodynamic mechanics of fish school. Previous physical models were passive flapping 
filaments that shed different vortices from real fish EDEl, or, rigid foils that have no 
body shape or body wave Ealia. Using a high-fidelity fish-like robot to uncover the 
hydrodynamic mechanism of the fish school has never been explored. 

In previous studies, on the energy harvest of hsh school and energy expenditure was 
mainly concerned in relation to lattices ll25]| . Among the formation hypotheses, a prevalent 
example is the diamond lattice first proposed by Weihs ifT^ . Recent studies, based on 
computational model, predicted that group configurations including diamond, rectangular, 
phalanx and line are all hydrodynamically advantageous ns El. Further studies on the 
real fish school demonstrated energy savings regardless of individuals’ spatial positions [D . 
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Whether spatial location is a crucial factor for energy expenditure when fish are swimming 
in groups is still controversial. Other factors related to energy consumption remain largely 
unexplored. 

Here, we consider a school of two fish-like robots swimming in a laminar flow. The 
high-fidelity of our robotic fish resides in its body shape, bio-inspired locomotion and sim¬ 
ilar Reverse Karman Vortexes (RKV) shedding by the caudal fin. The flow first encounters 
the rigid head and then the flexible body. Thus, a realistic body shape is critical. The loco¬ 
motion should be similar to the real fish to generate a similar flow field. Moreover, since 
the energy transfer between fish in school is mostly based on the vortices near the caudal 
fin, thus a similar wake is also significant. 

In this study, we test and compare the power consumed by two robot flapping fish fixed 
in a flow tunnel at different relative positions and flapping phase differences. To make fish 
body undulate in the flow field, the robot should consume power to overcome the resistive 
force of water. When robots condense into a group state, water-mediated interactions may 
lead to modifications of energy cost to keep robots undulating in the same body wave. 
Through series of experiments, we And that (a) the trailing fish will benefit energy saving 
when swimming in some conditions, and (b) when the schooling robotic fish swim in the 
similar spatial arrays of the real fish, the energy cost of the trailing fish mostly depends on 
the phase difference. 


2. Material and methods 

2.1. The high-fidelity robotic fish. The robotic fish used in our experiments is designed 
according to a kind of carangiform fish, which propel themselves by the undulation of the 
rear body and caudal fin. In order to improve the fidelity of the flow-mediated interactions 
between schooling fish, we enable our robotic fish to (a) have a controllable flexible body, 
(b) swim with biomimetic locomotion and (c) shed a similar flow field after the tail. 

The first step is to endow the robotic fish with a fish-like body. Based on the analyses 
of biologists, the flexible body has two almost mutually exclusive mechanical properties 
for swimming: flexible in lateral directions and incompressible in the longitudinal direc¬ 
tion ll26ll . In order to meet these requirements, we construct our robotic fish with three 
servomotors to generate the lateral rhythmic locomotion. The connections between them 
are rigid aluminium frames, in longitudinal director. The caudal fin is made with stiffness 
decreasing from the peduncle to the tail end, resembling to the real one. 

Most carangiform fish swim forward by passing a body wave down their bodies, a mode 
known as undulatory propulsion. One of the widely used body wave function is proposed 
by Lighthill l27l. Here we modified the function as Il2^ . 

(1) y = (ci(a — Ao) -|-C2 (a — Ao)^)sin(foc-|-2;r/f) 

where x is the displacement along the boy axis and y is the lateral body displacement (see 
Figure. [T](b)), Cl and C 2 are linear and quadratic wave amplitude envelope, k denotes the 
number of whole body wave length, / is the frequency of body wave, xq is defined as 
fixed point, which has no oscillation in y axis. The pattern of body undulation is shown 
in Figure, [^c). From the body wave function, we find that each segment of the fish body 
has its own oscillating amplitude as well as coupled phase differences with its adjacent 
segments. Guided by the body wave undulation, we implemented a similar undulation 
mode in our robotic fish by a bio-inspired controller. Central Pattern Generator (CPG) 
controller siMiini. The advantages of using CPG controller are that both the phase 
difference among the joints and between the fish can be well controlled on line. 
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Figure 1. An introduction of the robotic fish, including (a) three- 
dimensional mechanism consisting of rigid head, flexible body and cau¬ 
dal fin, (b) a snapshot of robotic fish swimming in water in the top 
view, (c) a cluster of body undulation derived from real carangiform fish, 
where BL represents body length, (d) a snapshot of vortices shedding by 
the robotic fish while swimming forward, (e) a typical schematic of the 
vortices, (online version in color.) 

Apart from a shape-similar body and bimimetic locomotion, the flow field shedding 
by the robotic fish is also critical for the interactions between “schooling” robotic fish. 
If we can generate a similar flow field around the robotic fish to the real one, the energy 
consumption of real fish school can be well estimated by our robotic fish. We compared 
the vortex shedding by real fish and our robotic fish to show the similarity between 
them EO). From the snapshot of the water flow of robotic fish in Figure. [T] (d), we can 
And vortices shed to the left-hand side of the robot rotate clockwise, and vortices shed to 
the right-hand side rotate counterclockwise (see Figure, [^(e)). This is quite similar to the 
Reverse Karman Vortex Street (RKVS), by which fish generates thrusts. 

2.2. Experimental Apparatus. To experimentally model the hydrodynamic interactions 
among fish school, we insert a “schooling” robotic fish into a laminar water tunnel, as 
shown in Figure. |^(a). The flow velocity of the water tunnel is set as f/ = 20 cm/s, 
which corresponds to the free-swimming speed of the robotic fish. The Reynolds number 
is 9 X 10^, which is similar to the natural situation. In this article, we study how the 
Lateral Spaces (LS), Streamwise Spaces (SS) and Phase Differences (PD) affect the power 
consumption for the tailing fish (see Figure. (b)). Among the three factors. Lateral 
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Lateral Spaces 


Flow 



Figure 2. An introduction of the test apparatus and the factors, (a) Lat¬ 
eral view of the experimental set-up. Two robotic hsh are fixed in a lam¬ 
inar flow with a free surface. A camera captures the interactions among 
the two robotic fish and the vortices, (b) Top view of fish swimming in 
group with Lateral Space (LS), Streamwise Space (SS) and anti-phase 
(Phase Difference (PD) equals n). (online version in color.) 

Spaces (LS) and Streamwise Spaces (SS) are controlled by group of step motors; and phase 
difference (PD) is controlled by the CPG controller ll30]l . Referring to the spatial lattices 
of real hsh school El and considering the limitation of the width of the How tunnel, we 
arrange the lateral space (LS) ranging from 0.27 to 0.33 BL, the streamwise space (SS) 
from 0.35 to 0.45 BL and the phase difference from 0 to 2n. 

As robotic Hsh is powered by electricity, a measurement and comparison of power con¬ 
sumptions of fish swimming alone and in a formation directly reflect the effects of swim¬ 
ming in groups. Robotic fish is powered by a stabilized volt supply (Matrix MPS-3005L- 
3). The electrical currents are measured by National Instrument Current Acquisition (NI 
9227) with a sample rate of 5000 Hz. The electrical current data acquisition lasts 10 s once 
after robotic fish swimming stably and repeats 5 times for each spatial lattice and phase 
difference. 

2.3. Evaluation of power consumption. Generally, the electrical energy input to the 
robotic fish will change to mechanical energy and other energy forms such as heat. For the 
swimming robot, the useful energy is the mechanical energy which is used to overcomes 
the reactive force of the fluid. We evaluate this part of energy consumption by subtracting 
the energy consumption when the fish is swimming in air from the energy consumption 
when fish swim in water 1331 iMl- Thus, the power consumption of fish swimming alone 
in water is evaluated as P„a — Pa, where P„a indicates the total power input when the fish 
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swims in water alone, and Pa represents the total power consumption while the hsh swims 
in air. Similarly, we calculate the power consumption for a robotic hsh swimming in school 
by Pws — Pa, where P„s is the total power consumption when the hsh swims in water in 
school. To compare the power consumption between hsh swimming alone and in groups, 
we dehne a power coefficient number rj. 


( 2 ) 



Pg) - {Pws 
Pwa - Pg 


Pa) 



A positive value of rj indicates a power reduction when hsh swimming in school. On the 
contrary, a negative value of J] means hsh in school will use more energy. Notably, all the 
values are time-averaged as the power consumption varies in one oscillation period. 

For each relative position and phase difference between robotic hsh, the tests include 
power acquisitions of robotic hsh swimming in air, in water alone and in a school. Uncer¬ 
tainties during the test quantities are dehned as the standard error of the mean a = a jy/N, 
for each situation consisting of N =5 measurements. And a is the sample standard devia¬ 
tion. 


3. Results 

3.1. Energy cost depends more on phase difference (PD). We hrst test the performance 
characteristics of a lone robotic hsh both in air and water. The power consumptions in four 
different conditions are shown in Figure, [^(a) (test in air). Figure, [^(b) (test in water). 
Figure, [^(c) (test in school with OAtz PD), Figure. |^(d) (test in school with 1.27r PD). 
Since the power consumption of robot is periodic at twice the frequency of body wave, the 
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Figure 3. Data analysis of power consumptions of the robotic hsh. 

The left column shows the time-domain analysis. And the right column 
shows the frequency-domain analysis. From the top to the bottom, data 
in each row are detected in different swimming conditions, (a) and (b), 
hsh swims alone in the air. (c) and (d) hsh swims alone in the water. 

(e) and (f) hsh swims in school with 0.4;r PD. (g) and (h) hsh swims in 
school with 1.2;r PD. 

power cost is evaluated by the average value. The average power cost of hsh swimming 
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in water P„a = 2.71 W is higher than that in air Pa = 1.44 W. Thus when fish swim alone 
in the water, the power cost to overcome the resistive force of water is P„a — fa = 1 -27 W. 
While fish swim in school with 0.4;r and 1.2;r, the power costs to overcome the resistive 
force of water P„s — Pa are respectively 1.32 and 1.21. Further, the efficiency coefficient 
for the 0.4;r PD is —3.94%, thus power increase. The efficiency coefficient for the 1.2;r 
PD is 4.72%, thus power reduction. 

Note that the noise signal is large in the power values acquired from the experiments. 
To verify the effectiveness of the data, we apply the Fourier transform to the data. As 
fish undulates the body periodically and the power consumption mainly depends on the 
reactive force of the fluid, the power cost value must also be periodic (as the filtered data 
shown in Figure, [^(a), (c), (e), (g)). Since the undulation of fish body is symmetrical, the 
frequency of the power 1.7 Plz is twice that of the body wave 0.85 Hz. Therefore, we prove 
the effectiveness of the data by this property, as shown in Figure. |^(b), (d), (f), (h). 



-4% I- e -^-1-^-. 

0 0.47r O.Stt 1.277 I.Ctt 277 

Phase Difference [Rad] 



0.28 0.3 0.32 0.28 0.3 0.32 

Lateral Space [BL] 


Figure 4. Phase difference is critical for fish swimming in groups to 
extract energy from the aquatic environment, (a) At the fixed Stream- 
wise Space (SS) 0.4 BL and each Lateral Space (LS), Power coefficient 
periodically varies across the whole phase. The shadow area shows the 
stand deviation for 5 trials at each point, (b) The power coefficient versus 
the spatial spaces at Phase Difference (PD) of 0.4 n. Most of the coef¬ 
ficients are negative denoting power increase, (c) The power coefficient 
versus the spatial spaces at Phase Difference (PD) of 1.2 n. Most of the 
coefficients are positive denoting power reduction, (online version in 
color.) 
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We then insert a second robotic fish into the flow with a fixed SS of 0.4 BL and series 
of LS. The power coefficient is calculated at each spatial lattice as well as phase difference 
for 5 trials, as shown in Figure. |^(a). The estimated error are illustrated in shadow. The 
variation of power coefficient as a function of PD is quite similar to a sine curve and are 
little related to spatial factors. Fish school swimming with PD around 0.4 n will consume 
more energy, in schools while around 1.2 ;r schooling will lead energy saving. Our data 
show that school fish energy saving is more related to phase differences than spatial lattice. 

To prove the hypothesis, we compare the power variation by fixing the phase difference 
0.47r Figure. |^(b), and 1.27r Figure. |^(c). In all the areas, fish swimming with phase 
difference of 0.47r costs more energy while phase difference of 1 .2n consumes less energy, 
indicating that the power consumption is more related to the phase difference than the 
spatial lattices. 



Phase Difference [Rad] 


Figure 5 . Two dimensional map of power coefficient for the follower 
swimmer in the (Phase Difference (PD), Streamwise Space (SS)) phase 
space. Lateral Space (LS) is discretized into (a) 0.27, (b) 0.29, (c) 0.31, 

(d) 0.33 (BL). The power reduction and increase are definitely divided 
by Phase Difference around 0.8 n. (online version in color.) 

To further study how the spatial and temporal factors affect the power consumption of 
group robotic fish. We carry out all the situations within the ranges, totaling 2600 trials. 
Figure. shows the power coefficient as a function of SS and PD with four different 
LSs. In contrast to swimming alone, the trailing fish may enjoy energy saving or suffer 
extra energy cost when swimming in different arrays and Phase Differences (PDs). For 
almost all the locations around the leader, the follower will benefit of saving energy when 
swimming with a phase difference ranging of 0.8;r ~L4;r. The larger the SS, the less 
energy the follower saves. However, we find that, at larger Streamwise Space (SS), the 
trailing fish need larger Phase Difference (PD) to get a maximum power reduction. This 
indicates that the spatial array also play a role in the power cost of the trailing fish but with 
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less impact than the Phase Difference (PD). Among all the trails we did here, the maximum 
power reduction for the follower is 6.9%, and the power increase is 5.6%. 

3.2. Hydrodynamic mechanism of energy saving in fish school. To understand why 
the follower in school saves energy and why the power consumption is strongly related 
to Phase Difference (PD), we carry out a techchnique of flow visualization to explore the 
hydrodynamics of fish swimming in school llTSl . We carry out the experiments by using 
laser-induced fluorescence (LIF) technique Il3^ . Since we would like to only show the 
interactions between the body wave and the shedding vortices, we put some fluorescent 
species on the fish body and inject the fluorescence dye from the tail. The fluorescent 
emission will be captured by the high-speed camera from the bottom of the tunnel at a fre¬ 
quency of 125 Hz (see Figure]^. As a result, only the parts with fluorescent species will be 
high lighted; other parts are black (Movies S4, S5). In order to illustrate the hydrodynamic 
mechanism of fish school clearly, we show the relationship between the vortices and the 
body wave in a schematic form in Figure. The left column (Figure. |^a), (c)) shows why 
fish in school saves energy in one period. And the right column (Figure, [^b), (d)) denotes 
the mechanism of why fish in school costs more energy in one period. The vortex dipole is 
shed by the leader and have an effect on the follower. The position of the robotic fish and 
vortices are all based on the Flow Visualization. 



Figure 6. Bottom view of the interactions between the vortices and 
undulation body of robots, (a), (c) The follower saves energy by flapping 
in the same direction of the vortex diploe induced by the leader, (b), (d) 

The follower costs more energy when flapping in the opposite direction 
of the vortex diploe induced by the leader (See Ref. ll^ for further 
explanation), (online version in color.) 

When the direction of the induced vortex dipole is the same as that of the flapping tail, 
the follower will benefits energy saving (see Figure. |^a), (c)). On the contrary, when the 
flapping direction is opposite to the direction of the induced vortex dipole, the follower 
will consume more energy (see Figure. |§b), (d)). Since PD describes the body wave of 
the follower leads or lags that of the leader, the relationship between the flapping direction 
and vortex-dipole direction is also determined by PD. When the schooling fish swim near 
(around 0.5 IBL), the energy consumption of the follower is greatly controlled by PD 
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instead of space array. For all the PDs, PD around 1.27r leads maximum energy saving for 
the follower swimming in school, and PD around QAn make the follower consume more 
energy. 

Further, one can dehnitely hnd that, with a phase difference of 1 .In, the follower swims 
in the same direction of the flow field induced by the vortex dipole, thereby saving energy. 
While the follower with a Phase Difference of QAn swims in the opposite direction, and 
thus it costs more energy. As the flow held induced by the vortex dipole is determined by 
the leader’s Happing tail. Therefore, the energy cost of the trailing Hsh essentially depends 
on the Phase Difference (PD), which proves our hypothesis. 


4. Discussion 

In this article, we applied a high-Hdelity robotic Hsh endowed with a controllable, flexi¬ 
ble and active body to investigate the hydrodynamic mechanism of power consumption in 
a Hsh school. Spatial and temporal factors were both considered in the robotic Hsh school. 
The power consumption of the trailing Hsh is tested and compared when Hsh are swim¬ 
ming alone and in schools. Moreover, we consider the spatial arrays that appeared in the 
real Hsh school Q. And under this spatial condition, we vary the Phase Difference (PD) 
between the two robotic Hsh swimming together. After series of experiments, we And that, 
under the same spatial condition as in real Hsh school, the trailing robotic Hsh saves energy 
depending more on the PD rather than the array spacing. 

Compared with the previous apparatuses, the robotic Hsh consists of flsh-like body, 
bio-inspired locomotion and similar Reverse Karman Vortex Street (RKVS). A popular 
apparatus that was used to uncover the mechanism of energy saving in Hsh school is the 
flexible filament in a flowing soap Him ETIImI . As the filament is flexible, the interactions 
between the body and fluid vortices can be well illustrated i32l. However, because of 
the passive body, the vortices shedding by the filaments are Karman Vortices instead of 
Reverse Karman Vortices. Using active flapping foils, B. M. Boschitsch and P. A. Dewey 
respectively studied the hydrofoils in a in-line configuration ll2^ and side-by-side m 
configuration. Although the rigid foils are active and controllable. They are quite different 
from the flexible tails, as real fish tail has a increasing elasticity towards the end. Moreover, 
a flapping foil can not represent the hydrodynamics of fish swimming in water, as the flow 
first encounters the rigid head and then interacts with the flexible body. The robotic fish 
we applied here combines both the advantages. And thus, using robotic fish to mimic the 
fish school to explore the mechanism of energy saving is feasible and more realistic. 

Note that the maximum energy saving for the follower is only 6.9%, smaller than pre¬ 
vious studies Gniminiiia. One explanation of this is that the real school is more than 
two fish. In this article, we only considered the basic schooling of two fish. As one leader 
will make the follower save 6.9%, two leaders will make the follower save about 13.8%. 
Moreover, the frequency of our robotic fish is around 3 times lower than the real fish, there¬ 
fore the energy of the vortex field shed by our robotic fish is smaller than that shed by the 
real fish. This may also lead our data show that the robotic fish only gets a small energy 
saving. But in the filed of engineering, an energy saving of 6.9% is great for the manmade 
machines. Our finding in the study may guide the formation and locomotion control for 
group of bioinspired underwater automatic vehicles. 

The work presented here is basis and needs more work. The power costs of the leader 
will be further explored. Nevertheless, this work first implements artificial fish school, 
and evaluated the mechanism of fish swimming in groups in the aspect of hydrodynamics. 


10 


L. LI, L. JIA, AND G. XIE 


Further the connections of biology and physics also lie in other aspects such as group 
behaviour verifications by using bio-inspired high-fidelity robots. 
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